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A theoretical model for the reentrant charge-order transition in the manganites is examined. This 
transition is studied with a purely electronic model for the e g Mn electrons: the extended Hubbard 
model. The electron-phonon coupling results in a large nearest-neighbor repulsion between e g elec- 
trons. Using a finite-temperature Lanczos technique, the model is diagonalized on a f 6-site periodic 
cluster to calculate the temperature-dependent phase boundary between the charge-ordered and 
homogeneous phases. A reentrant transition is found. The results are discussed with respect to the 
specific topology of the f6-site cluster. 



The manganites have a very rich phase diagram that 
includes ferromagnetic, antiferromagnetic, and charge- 
ordered phases Various theoretical models have 
been used to explain different aspects of this phase dia- 
gram {§-§. 

In its simplest incarnation, the charge-ordered (CO) 
phase occurs at hole doping x = 1/2 with equal amounts 
of Mn 3+ and Mn 4+ ordered in real space in a checker- 
board pattern. The oxygens relax away from the Mn 3+ 
ions and towards the Mn 4+ ions, thus providing a re- 
pulsive potential between Mn 3+ ions (or equivalently be- 
tween Mn 4+ ions). The potential energy gain exceeds the 
kinetic energy loss due to the formation of this insulating 
state M. 

When observed, the CO is generally the lowest 
temperature phase, but recently the CO phase has 
been seen to melt with decreasing temperature in 
Pro.65(Cao.7Sro. 3 )o.35Mn0 3 @ and LaSr 2 Mn 2 7 @. 
The lowest temperature phase is metallic, and the CO in- 
sulator is only observed at intermediate temperatures. A 
reentrant transition has been obtained theoretically using 
extended Hubbard models both with electron-phonon in- 
teractions JlT| and without electron-phonon interactions 

& 

In this paper, we study the charge-order transition in 
the extended Hubbard model (without electron-phonon 
interactions) on the two-dimensional square lattice. Pre- 
vious work [jl2j solved this model in infinite spatial di- 
mensions, resulting in hnite entropy (due to the spins) 
at T = in the CO phase, so a reentrant transition was 
guaranteed to be found. In the infinite two-dimensional 
square lattice, the spins will order into a Neel state with 
zero entropy at T = 0. 

The Hamiltonian is given by 

H = t i c \a c o<y + n - c + U X! n *T n U + Vy^^jTij, 

(ij)o- i (ij) 

(1) 

where c^ CT (c$ CT ) creates (annihilates) an electron with spin 
cr on site i, is the number operator with spin a on site 
i, and m = rii\ + n^. The hopping amplitude is t, (ij) 



enumerates nearest neighbor sites on the two-dimensional 
square lattice, U is onsite repulsion, and V is the nearest- 
neighbor repulsion. The non-interacting bandwidth on 
the two-dimensional square lattice is given by W = 8|i|; 
we set U — W and vary V at quarter filling (one electron 
for every two sites). For small V, we expect the ground 
state will be a homogeneous Fermi liquid. For large V, 
the electrons will crystallize in a checkerboard pattern to 
avoid occupying neighboring sites. 

We solve the Hamiltonian ([!]) on a 4 x 4 cluster using a 
recently developed hnite-temperature Lanczos technique 
0,0. We choose periodic boundary conditions result- 
ing in a closed Fermionic shell in the non-interacting limit 
Jl5| . In each symmetry sector, we perform Nl = 4000 
Lanczos steps. We eliminate spurious eigenvalues [ jl6| , 
leaving more than 2000 real eigenvalues in each sector, 
but not all of these will be converged. The extreme 
eigenvalues (lowest and highest) converge hrst [O] . From 
these eigenvalues, we compute the susceptibility 



X ~NdV' 



(2) 



where N — 16 is the number of sites. The susceptibility 
X is equivalent to the nearest-neighbor pair correlation 
function by the Kubo formula. In the homogeneous phase 
X is large, but in the CO phase where there is only a small 
probability of finding electrons on neighboring sites, x is 
small. 

The susceptibility (^) for five different temperatures is 
plotted in Fig. ^. For all temperatures shown, x 1S large 
at small V in the homogeneous phase, decreases contin- 
uously with increasing V, and flattens out at a small 
value in the CO phase at large V. To determine the 
boundary between the homogeneous and CO phases, we 
pick a critical value of Xc = 0.1 for the phase boundary. 
We choose this value to coincide approximately with the 
maximum of d 2 x/dV 2 . A different value of Xc results in 
a phase boundary that is qualitatively the same: for all 
reasonable values of Xa the phase boundary first shifts 
to smaller V as the temperature is raised from zero and 
then moves to larger V in the high-temperature regime, 
characteristic of a reentrant transition. 
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FIG. 1. Temperature dependence of the susceptibility as a 
function of the nearest-neighbor repulsion V (in units of the 
bandwidth W = 8|£|) for the 4x4 periodic cluster at U = W. 
The susceptibility varies continuously even at T = (solid 
curve). The susceptibility is large in the homogeneous phase 
(small V) and small in the CO phase (large V). With in- 
creasing temperature, the susceptibility first shifts to smaller 
V then to larger V, indicating a reentrant transition. 
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FIG. 2. Phase diagram for the model with Xc = 0.1. 
There is clear reentrant behavior around V/W « 0.35. The 
four curves correspond to different numbers of Lanczos steps. 
Runs with more Lanczos steps result in more eigenvalues, and 
thus are accurate to higher temperatures. The Nl = 4000 
run is clearly accurate well into the high-temperature homo- 
geneous phase in the reentrant region. 



The phase diagram calculated from the point where 
x(V,T) = Xc = 0.1 is shown in Fig. g. Reentrant behav- 
ior is seen for 0.31 < V/W < 0.37. For V in this range, 
the ground state is homogeneous, but a CO phase exists 
at intermediate temperatures. 

The 4x4 system has no entropy in cither the homoge- 
neous or the CO phases, so the phase boundary in Fig. 
^ has infinite slope in the T — » limit. However, the 
CO phase has more low-lying excitations than the homo- 
geneous phase, so in the reentrant region the CO phase 
becomes favored at intermediate temperatures. 

The greater number of low-lying excitations in the CO 
phase may be at least partially due to the topology of the 
4x4 periodic cluster used. The homogeneous phase will 
have the usual Fermi-liquid excitations. In the CO phase, 
the electrons order in a checker-board pattern. The spins 
interact via a fourth order superexchange process 0,|l8| ■ 
A given spin can interact both with its diagonal neigh- 
bor and with its neighbor two steps away either horizon- 
tally or vertically. On the infinite lattice, the interaction 
with the diagonal neighbor is 4 times as strong as the in- 
teraction with the horizontal/ vertical neighbor because 
there are 4 times as many diagonal superexchange paths. 
Thus the spin system is equivalent to the extended anti- 
ferromagnetic Heisenberg model on a square lattice with 
J i = 4J 2 , where J\ is the nearest- neighbor interaction 
and Ji is the next-nearest-neighbor interaction. Exten- 
sive numerical calculations have shown that this Heisen- 



berg model forms a Neel state at T = 19 2(|. 

Because two steps to the left is the same as two steps 
to the right (and the same for up and down), the 4x4 pe- 



riodic cluster is equivalent to a hypercubic 2 4 cluster [^l| . 
Thus the diagonal superexchange coupling is the same as 
the horizontal/ vertical coupling. So the CO Neel state on 
this particular cluster is strongly frustrated and may not 
form. This could contribute to the greater density of low 
lying states in the CO phase than in the homogeneous 
phase. 

Smaller periodic clusters have even more severe finite- 
size problems. We decided to work with the 4x4 clus- 
ter because it is the largest system with a manageable 
Hilbert space. Even with all symmetries applied, the 
largest Hamiltonian matrices were larger than 2.0 x 10 s 
(the matrices are very sparse, with only about 30 non- 
zero elements in each row or column). The Hilbert space 
diverges exponentially with the number of sites in the 
cluster. The next largest cluster with an even number 
of electrons has 20 sites and has irreducible sectors with 
Hilbert spaces larger than 2.4 x 10 8 , more than 1000 times 
larger than those of the 16-site system. 

In summary, we have examined the reentrant charge 
order transition in a simple model for the manganites, 
the extended Hubbard model in two dimensions. We 
computed the low-lying eigenvalues of the model on a 
16-site periodic cluster, and used the susceptibility with 
respect to the nearest-neighbor repulsion to compute the 
low-temperature phase boundary between the homoge- 
neous and charge-ordered phases. Like previous results 
on this model in infinite dimensions, we find a parameter 
region where the model shows reentrant behavior. The 
reentrant region may be partially an artifact of the small 
cluster used, which frustrates the formation of Neel spin 
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order in the charge ordered state. 
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